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Abstract. We present a new, short proof of the increased regularity obtained 
by solutions to uniformly parabolic partial differential equations. Though this 
setting is fairly introductory, our new method of proof, which uses a priori 
estimates, can be extended to prove analogous results for problems with time- 
dependent coefficients, transport equations, and nonlinear equations even when 
other tools, such as semigroup methods or the use of explicit fundamental 
solutions, are unavailable. 


1. Introduction. It is well-known that solutions of uniformly parabolic partial 
differential equations possess a smoothing property. That is, beginning with ini¬ 
tial data which may fail to be even weakly differentiable, the solution becomes 
extremely smooth, gaining spatial derivatives at any time t > 0. Though this prop¬ 
erty is well-established, such a result is often excluded from many standard texts 
in PDEs [2,3,5- 7,11]. Of course, these works contain theorems demonstrating the 
regularity of solutions, but the same degree of regularity is assumed for the initial 
data. The notable exception is [1, Thm 10.1] in which a gain of regularity theorem 
is proved, specifically for the heat equation with initial data in L 2 (fl) using semi¬ 
group methods. Even more concentrated works on the subject of parabolic PDE 
[4,8] do not contain such results regarding increased regularity of solutions to these 
equations. 

In the current paper, we will present a few results highlighting the increased 
regularity that solutions of these equations possess. In particular, our method of 
proof is both new and brief, and relies on a priori estimates. Hence, when traditional 
tools like semigroup methods or explicit fundamental solutions cannot be utilized, 
the new approach contained within may still be effective. 
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2. Main Results. Let n £ N be given. We consider the Cauchy problem 

j dtu — V • ( D(x)S7u ) = f(x), x £ K", t > 0 
u(0, x) = uo(x), 


where D , the diffusion matrix, and /, a forcing function, are both given. Equations 
like (1) arise within countless applications as diffusion is of fundamental importance 
to physics, chemistry, and biology, especially for problems in thermodynamics, neu¬ 
roscience, cell biology, and chemical kinetics. As we are interested in displaying the 
utility of our method of proof, we wish to keep the framework of the current prob¬ 
lem relatively straightforward. Thus, we will assume throughout that the diffusion 
matrix D = D{x) satisfies the uniform ellipticity condition 


w ■ D(x)w > 9\w\ 2 


( 2 ) 


for some 9 > 0 and all x,w £ R". We note that under suitable conditions on the 
spatial decay of u, our method may also be altered to allow for diffusion coefficients 
that are not uniformly elliptic (see [10]). Additionally, we will impose different 
regularity assumptions on D and / to arrive at different conclusions regarding the 
regularity of the solution u. 

Throughout the paper we will only assume that the initial data uq is square inte¬ 
grate. Hence, even though uq(x) may fail to possess even a single weak derivative, 
we will show that u(t,x ) gains spatial derivatives in L 2 on (0, oo) x R n . Hence, by 
the Sobolev Embedding Theorem, solutions may be classically differentiable in x 
assuming enough regularity of the coefficients. In addition, we will show that u is 
continuous in time at any instant after the initial time t = 0. Though the setting 
(1) is fairly introductory and the assumptions on D and / are somewhat strong, the 
new method of proof can be adapted to extend the results to problems with time- 
dependent terms, transport equations, systems of parabolic PDEs, different spatial 
settings such as a bounded domain or manifold, and nonlinear equations, includ¬ 
ing nonlinear Fokker-Planck equations and nonlinear transport problems arising in 
Kinetic Theory [9,10,12], 

For the proofs, we will rely on a priori estimation and the standard Galerkin 
approximation to obtain regularity of the approximating sequence and then pass to 
the limit in order to obtain increased regularity of the solution. Hence, we focus 
on deriving the appropriate estimates as the remaining machinery is standard (cf. 
[2,5]). In what follows, C > 0 will represent a constant that may change from line 
to line, and for derivatives we will use the notation 



to sum over all multi-indices of order k £ N. When necessary, we will specify 
parameters on which constants may depend by using a subscript (e.g., Ct)- 
Our first result establishes the main idea for low regularity of D and /. 


Theorem 2.1 (Lower-order Regularity). Assume f £ 1L 1 ( R"), D £ kF 1,00 (R"; R"^" 
and uq £ L 2 (R"). Then, for any T > 0 and t £ (0,T], any solution of (1) satisfies 


SUp \\u(t)f 2 < C T (\\U 0 \\1 + ||/||hi) and l|V x u(f )||2 < ^(|K||I + ll/llffi)- 
<t<T t 
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Proof. We first prove two standard estimates of (1). First, we multiply by u, inte¬ 
grate the equation in x, integrate by parts and use Cauchy’s Inequality to find 

|4ll«(*)|||+ / v x u{t) ■ DV X u{t)dx < i (II/II 2 + \\u(t)\\l) . 

J R n 


Then, using (2) and the regularity assumption on /, we find 

“Nt)ll! < c (||/|| 2 ffl + ||«(t)||i) - 0||v x «(i)||!. (3) 

Next, we take any first-order derivative with respect to x (denoted by d x ) of the 
equation, multiply by d x it, and integrate to obtain 

\^:\\d x u{t)\\\+f V x d x u(t)-DV x d x u(t)dx + [ V x d x u-d x D\7 x u dx = f d x fd x udx. 

2 dt Jr" Jr" Jr™ 

Thus, using Cauchy’s inequality with the ellipticity and regularity assumptions, we 
find for any e > 0 

\-jl\\dxu{t)\\l < - l V x d x u(t) • DV x d x u(t)dx - f V x d x u ■ d x DS7 x u dx + [ d x fd x u 
* dt Jw 1 ir J R n 

< -9\\V x d x u{t)\\l + \\D\\ w i,<*, (e\\V x d x u{t)\\l + i||V x u(f)||^ 

+ \ (\\9xf\\l + \\d x u{t)\\l) . 

Choosing e = 9{2\\D\\ w i, ao )~ 1 and summing over all first-order spatial derivatives, 
we finally arrive at the estimate 

\j t \\V x u{t)Wl < C (ll/ll^i + ||V x «(t)||l) - \\\Vlu{t)Wl (4) 

Now, we utilize a linear expansion in t to prove the theorem. Let T > 0 be given. 
Consider t € (0,T] and define 


Ml (t) = \\u(t)\\ 2 2 + -\\V x u 


We differentiate this quantity, and use the estimates (3) and (4) to find 


M[(t) = -\\u(t)\\ 


9 „ , , ll0 9t d „ 

^ Mt )\\ 2 2 + —^-l|V x u( 


< 


(ii/n 2 ffl + nmi) - m\vMt)\\i +^iiv xU (t)ii 


9t r 

”2 r 


< C T 


2 C 


H 1 


2 

H 1 


Mt) 111) -s\\v 2 x u{ 


+ Afi 


A straightforward application of Gronwall’s inequality (cf. [2]) then implies 

Mi(t) < CAM.i 0) + WfWU = C T (\\uo\\l + \\ff Hl ). 

Finally, the bound on M\ (t) yields 

IK*) 111 < C T , ||V x «(t)||| < ^ 

and the estimate holds on the interval (0,T]. As T > 0 is arbitrary, the result 
follows. 

□ 
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Next, we formulate the existence of weak solutions for our lower-order regularity 
setting. 

Definition 2.2. We say that u £ L 2 {\ 0, T]; H 1 (W 1 )), with dtu £ L 2 ([ 0, T]; 
is a weak solution of (1) if 

f (d t u, v) + (D(x)\7u, Vu) = (/, v) 
u(0,x) = Uq(x), Itl" 

for every v £ i? 1 (K") and t £ [0,T]. 

Theorem 2.3 (Existence and Uniqueness of Weak solutions). Given any uo £ 
L 2 (R n ) and f £ H 1 (M. n ) and T > 0 arbitrary, there exists a unique u £ C(( 0, T]; H 1 (R n ))n 
C([0,T];L 2 (K”)) and u' £ T 2 (([0,T ];that solves (5). 

Proof. We follow a standard Galerkin approach. Take {w)fc(x)}^U 0 to be an or¬ 
thonormal basis for L 2 with Wk £ H s for s > 0. Consider functions of the form 
u m (x, t) = Yk=o dk(t) w k{%) with d™{t) a smooth function of t. Then the equations 

{dtu m , w k ) + (D(x)X7u m , Vw k ) = (/, w k ) 

(u m (0,-),w k ) = (u 0 ,w k ) 

for k = 1, 2,... m reduce to a constant coefficient first order system of ODE’s for 
d™(t), and hence existence of approximate solutions is readily established. 

For these solutions, u m (t ), we may repeat the proof of our a priori estimates 
verbatim. Thus we can conclude that 

( Qt \ 

\\Um{t)\\l + TT \\V x U m (t)\\l ) < C T llnollz ■ 

From the proof, we also have the inequality 

j 0 \jt dt + j Q dt < C T (||/||2 + IKH2 ) dt. 

Using the above control of sup 0<t<T ||w. m (f)||i, we find that 

f ||Wm||tfl( R n) < Cr( H/H 2 + 11^0 II 2 ) ' 

J 0 

Finally, fix v £ with ||v|| ff i < 1 and consider 

(dtu m , v) = - { D(x)\7 x u m , V x i>) + (/, v). 

Using Cauchy-Schwartz and taking the supremum over v £ H 1 with ||u|| ff i < 1 we 
find ||d t Um(*)|| ff -i < ||/|| 2 + ||«o|| 2 ) and 

J \\ d tUm(t)\\ 2 H -i dt<C T (\\f\\l + \\uo\\ 2 2 y 

Thus, u m is a bounded sequence in L 2 ([0, T]; Ff 1 (K")) and d t u m is a bounded 
sequence in L 2 ([0, T]; fF _1 (R”)), so we may extract a subsequence rrij so that u mj —> 
u in L 2 ([0,T]; iF 1 (K")) and dtu mj —> d t u in L 2 ([0, T]; H~ 1 (W n )). 

Now fix an integer N and consider v(t) = Y^=o dk(t)wk{x), where dk(t) are fixed 
smooth functions. Then for mj > N, we have 
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Thus passing to the limit 



(d t u, v) + {D(x)Vu 1 S7v) dt 



(f,v) dt. 


Since v given above are dense in L 2 ([0, T]; i? 1 (R”)), the equality holds for any v 
in this space. Since u £ L 2 ([0, T];H 1 (W 1 )) and d t u £ L 2 ([0, T]; H~ 1 (W n )) we have 
that u £ C([0,T];L 2 (R n ) by [2, Thm. 3 §5.9.2], 

To see m(0) = uq, we take v £ C' 1 ([0, T]\H 1 (M. n )) with the property that v(T) = 
0, then we find that 


(«m,(0),v(0)) 



{u mj ,dtv) + {D(x)Vu mj ,v) dt 



(/, v) dt. 


Notice u mj { 0) = ( u o, w k) Wk —t uq in L 2 (R”) as nrij —> cso. Thus passing to 

the limit, we find that (u(0),u(0)} = (uo,u(0)) for u(0) arbitrary. Hence u( 0) = uq . 

To prove uniqueness, notice that for any two solutions u and u the difference 
u — u satisfies our equation with u 0 = 0 and / = 0. Thus our a priori estimate gives 
that sup 0<t<T ||u(f) — fi(t)H 2 < 0 and uniqueness follows immediately. 

Finally, to show that u £ C{{ 0, T]; H 1 (W n )) we consider w s (t) = u(t + s) — u(t). 
Then w s (t) satisfies our equation with / = 0 and w(0) = uq — u(s). From our a 
priori estimate 


Of 

KWH + j \\V x w e ml < C T (\\u 0 - u(s)\\l). 

From the fact that u £ C([0, T]; L 2 (R")), we have for t > 0 that lim s _u) \\u(t + s) — u(t)\\ 
0 and lim s _> 0 || W x u{t + s) — V x M(t)|| = 0, whence the result follows. □ 


Next, we extend the previous estimate to higher regularity assuming that D and 
/ possess additional weak derivatives. 


Lemma 2.4 (Higher-order Regularity). For every m £ N, if f £ D £ 

kF rn,00 (M”; M" XTl ), and uq £ L 2 (R"), then the previously derived solution of (1) 
satisfies 

\\ V x u \\l - -fir (KII 2 + ll/llff*) for k = 0,1,..., m. 

Proof. We will prove the result by induction on to. The base case (to = 1) fol¬ 
lows immediately from Lemma 2.1. Prior to the inductive step, we first prove a 
useful estimate for solutions of (1). For the estimate, assume D and / possess 
k £ N derivatives in L°° and L 2 , respectively. Take any /cth-order derivative with 
respect to x (denoted by <9“) of the equation, multiply by and integrate using 
integration by parts to obtain 


li\\ d X t )\\l + f ■ ]T ]T ( a ')d?D\7 x d2udx= f dfffdffudx 

zat JR" ,-_n \PJ J R" 


1 =° 101 =3 
0+7=a 


and thus 


I V x dfu ■ Yi Y (Z] d x D ^xd2u dx+ [ dfffdffu dx. 
JR" -_ n \PJ JR" 


1 =° 101=1 
0+7=a 
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Labeling the first term on the right side A, we use (2) and the regularity of D 
with Cauchy’s inequality (with e > 0) to find 


A = 


< 

< 



V x d x u ■ DV x d x ii dx — 






/3+7=a 


-9\\\7 x d:u{t)\\l + C\\D\\ 2 wKoo ^||V x S>(t)||2 + i||n(t)||^_ 1 ) 

- 0 -\\v x dymt+c\\um 2 Hk -i 


where we have chosen e = 8 (2C'||D||^ fci00 ) 1 in the third line. Inserting this into 
the above equality and using Cauchy’s inequality again we find 

~\\d>m 2 2 < -^l|V^>(i)||| + C\\u(t)\\ 2 H^ +l\\dSf\\l + 

Finally, summing over all first-order derivatives and using the regularity of / yields 
the estimate 


~||V*«(t)ll 2 < -|||V*+ 1 «(*)||1 + c ( W.f\\ 2 Hk + lKt)|&„). (6) 

Now, we prove the lemma utilizing this estimate for k = 0,1, m. Assume / £ 
H m (R n ) and D £ bF m ’ 00 (M"; Since this implies that / £ H m ~ 1 {W l ) and 

D £ W™- 1 ’ 00 (M. n ; M nxra ), we find that u £ C((0, oo); by the induction 

hypothesis. Let T > 0 be given. Consider t £ (0,T] and define 


M(t) = y 

fc=o 


(9t) k 
2 k k\ 


I Vtui 


2 

2’ 


We differentiate to find 


m nk-i.k—1 

M ® = ^ 2 k (k — 1 )! 


w^Mmi+Y 


fc=o 


(. 9t) k d 
2 k k\ dt 


\V k u(t)\\ 2 2 =:I + II. 


We use (6) for any k = 0,..., m and relabel the index of the sum so that 
( 0t) k 


fc=0 


11 ^ 

m (s\j.\k 

+ C E ^Fk\ {^^ Hk + 


m pk-\- lj.k 

= - 2 £L^||v^V 


k =0 


2 k+1 k\ 


k =0 


—2/ — " ~ , ||V™ +1 r 

2 m + 1 m! 11 


Will + jejj (ll/lla* + II“( () IIh.) 


k—0 


< 
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Notice the induction hypothesis gives \\u(t)\\ 2 H k-i < p&GI/Hffk-i + ILoll^)- Using 
this bound and the previous inequality within the estimate of M'(t), we find 


TTl / /)j-\k 

M '(t ) < c£^ r (ll/ll^ + K*)lll.) 

k =0 

< C T (\\ff Hm +it { ^ [Il|V*«(i)||l + HtiWII^-x] ) 

V fc=0 / 

< Ct (\\f\\ 2 H m + M(t) + ^2 ^fcTj-TfcZi (ll/llfffe-1 + WuoWl)) 

V fc=0 ' / 

< C T {\\f\\ 2 Hm + \\u 0 \\l + M(t)). 


Another straightforward application of Gronwall’s inequality then implies 

M(t) < C T (\\ff Hm + \\u 0 \\l + M{ 0)) < C T (\\ff Hm + Huolll). 

Finally, the bound on M(t) yields 

l|V™«(*)lll < pp 

which completes the inductive step and the proof of the lemma. □ 

Theorem 2.5 (Existence and Uniqueness of Weak solutions). Let m € N be given. 
For any uq £ L 2 (M. n ), f £ H m ( R”), and T > 0 arbitrary, there exists a unique 
u e C((0, T\; 7J m (R"))nC'([0, T]; L 2 (R n )) and v! £ T 2 (([0,T]; H~ 1 (M n )) that solves 
(5). 

Proof. This follows by a straightforward repetition of the proof of Theorem 2.3 with 
the obvious modifications. □ 

Of course, if the dimension n satisfies n < 2m — 1 this result implies classical 
differentiability of solutions and they satisfy the PDE in the classical sense. Finally, 
this result can be easily used to deduce infinite spatial differentiability of the solution 
assuming D and / satisfy the same condition. 

Theorem 2.6 (Infinite Differentiability). If f £ H°°(R n ), D £ VF 00,00 (IR n ; K raxn ), 
and uq £ L 2 (K”), then for any T > 0 arbitrary, any solution of (1) satisfies u £ 
C°°((0, T] x R"). 

Remark 1. On a bounded domain, it is enough to impose / £ C°°(R. n ) and 
D £ C'°°(R";R" xrl ) to arrive at the same result. 

Proof. The result follows immediately by applying Theorem 2.4 for each to £ N, 
noticing that 

d t u=V- (DWu) + f 

is continuous, and bootstrapping this property for higher-order time derivatives. □ 

Remark 2. Though we have chosen to demonstrate the method for equations with 
time-independent coefficients, the same results can be obtained for time-dependent 
diffusion coefficients D and sources / using the same proof, as long as these functions 
are sufficiently smooth in t. Additionally, similar arguments can be used to gain 
regularity of the solution in t, as well. 
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